Spectral properties of the one-dimensional two-channel Kondo lattice model 
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We have studied the energy spectrum of a one-dimensional Kondo lattice, where the localized mag- 
netic moments have SU(N) symmetry and two channels of conduction electrons are present. At half 
filling, the system is shown to exist in two phases: one dominated by RKKY-exchange interaction 
effects, and the other by Kondo screening. A quantum phase transition point separates these two 
regimes at temperature T = 0. The Kondo-dominated phase is shown to possess soft modes, with 
spectral gaps much smaller than the Kondo temperature. 
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In this paper we discuss a one-dimensional two-channel 
SU(N)-symmetric Kondo lattice, that is, a 1-D lattice 
where spin operators are represented by generators of 
the su(N) algebra and the conduction band has an ad- 
ditional two-fold degeneracy. We discuss the case where 
each conduction band is half-filled. At this point due to 
the additional particle-hole symmetry the total symme- 
try of the model is SU(N)xZ 2 xSU(2). 

The interest in this model is justified by the fact that 
the corresponding impurity model exhibits non- Fermi li- 
quid (NFL) behaviourcl, which alsa takes place for the 
D = oo two-channel Kondo latticea. The origin of this 
NFL behaviour in the single impurity problem is quite 
clear - it is caused by dispersionless modes (see Refs. [H] 
)■ 

Our study of the half-filled lattice case demonstrates 
that the system can exist in two phases: the RKKY 
dominated phase, where, its low-energy behaviour is de- 
termined by the RKKYq exchange interactions between 
local magnetic moments induced by the conduction elec- 
trons, and the Kondo-dominated or heavy fermion state, 
where Kondo screening of the magnetic impurities takes 
place. Both phases are separated by a quantum phase 
transition point at T = 0. This agrees with_the scenario 
for the Kondo lattice suggested by DoniacbO. The heavy 
fermion phase possesses soft modes; these modes have 
spectral gaps which are much smaller than the Kondo 
temperature (Tr). 

Most of the papers devoted to the Kondo lattice use the 
slave-boson approach and the large-N approximationlj, 
which is more suitable for the Kondo dominated (disor- 
dered) phase. The opposite situation is easier to study in 
one dimension, where the RKKY interaction has a sim- 
pler form. This has been done by one of the authors 
for the one-dimensional M-channel Kondo lattice at half- 
fillingu: when the spin sector has SU(2) symmetry. In 
this case both charge and orbital sectors acquire gaps, 
while for half-integer value of |M/2 — S\ the spin sector 
remains gapless and at low energies is described by the 



SUi(2) Wess-Zumino-Novikov-Witten model (WZNW)i. 
In this paper we generalize this analysis for SU(N) sym- 
metry of the spin sector: the existence of an additional 
parameter, N, opens a possibility to study a competition 
between the RKKY- and the Kondo-dominated regimes. 
We consider the following Hamiltonian: 

i;m=l,N;a=l,2 

+ J C tm, a d i,mdt tl Ci t l ia +H.C. (1) 

i;m,l=l,N;a=l,2 

Here i is the lattice site index, m{— 1, ...N) numbers com- 
ponents of the SU(N) spin, while a is the channel index 
of the conduction electrons. cf m a is the creation opera- 
tor for conduction electrons in a state characterized by 
quantum numbers i,m,a>, while df m corresponds to the 
local moment states. The Hamiltonian is supplemented 
by the constraint 

df, m di,m = q, (2) 

m=l,N 

where the number q represents the filling of the local 
orbitals. It remains finite in the limit N — ► oo. 

The calculations for the RKKY regime are similar to 
the ones for the SU(2)-symmetric case done inQ. The 
mean field gap in the charge sector is A ~ J, while in 
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the Kondo regime it is Tk ~ De NJp ^F) ; being D the 
conduction band bandwidth and p{cf) ~ D^ 1 the den- 
sity of states at the Fermi level. For large N, the latter 
gap is the largest one. 

The appearance of the gap leads to the change of the 
ground state energy; in the RKKY-phase we have 

Np(e F )A\ D Np(e F )J 2 D 
Erkky = m-r « In — (3) 
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and in the Kondo regime it is: 
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(4) 



The critical value of TV at which both ground state ener- 
gies become comparable is colossal for any realistic J: 



Jp{e F ) 



In-^l/Jpiep)} 



(5) 



but it does not matter since we are going to study uni- 
versal properties of the spectrum. 

Let us now discuss the properties of the system on 
both sides of the transition. Deep in the RKKY phase 
the adequate description of the spin sector in the con- 



non-linear 



tinuum limit is given by the u[N/2)m[N/2) 
sigma model (NLct) with the topological termil. In the 
limit TV = 2 this becomes the SU(2)/U(1) modeUwith 
the topological term obtained for this case in Ref . [0] . In 
the presence of the topological term the model is likely to 
be massless. The transmutation scale (that is the energy 
scale at which the model undergoes a crossover to the 
strong coupling regime) of this model is: 
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g ~ p(e F )Jy/ln(D/J). 



(6) 



On the other hand, in the Kondo limit the spin sector 
most certainly has a gap (we shall see that it is much 
smaller than TV). Since the NLer model is critical, its 
renormalization group beta functionta has a zero at some 
finite value of g = p(ep)JN (see Fig.l). If there is no 1st 
order phase transition between the two regimes, the beta 
function can be continued into the Kondo regime where 
the spectrum has a gap. For such behaviour to take place, 
one needs to have a beta function which does not cross 
the horizontal axis but just touches it. 



RKKY 



Kondo 




Fig. 1. Schematic plot of the beta function as a function 
of the bare coupling. 



This behaviour of the (3 function corresponds to the 
presence of a marginal operator at the critical point. 

We will now discuss the Kondo-dominated (large TV) 
regime in detail. The Lagrangian density can be written 
as follows: 



i ,j;m=l,N;a = l,2 

+ E d tm d r di 



i\m=l,N 
J 
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i:m,l=l.N:a=l,2 



(7) 
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In this expression A is a Lagrange multiplier introduced 
to enforce the local constraint on occupation number of 
the localized orbitals. 

We next decouple the interaction term employing the 
Hubbard-Stratonovich transformation: 
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C- ™ ^di m d^ ^C-i^l^a. 
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E [Vi,act m , a di, m + H.c], (8) 

i;m=l,N;a=l,2 



where V itCt =\ V | z i<a and J2 a =i,2 z ta z h» = 1 - Wc 
now apply the following further transformation for the 
conduction band orbitals: 

which can be cast also in the following matricial form: 
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(10) 

Here g is a matrix belonging to SU(2). 

One can imagine two simple-minded mean field confi- 
gurations now: 

a) the uniform one, in which the g matrices on different 
sites are equal; 

b) the case Z2i+\. a — £ a ,p b > which in matrix form 
would read: 



ff2i+l — — 32i 
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In the first case, in mean field approximation only one 
conduction band channel (x-chain) hybridizes with the 
local orbitals. Meanwhile, in case (b) both chains are 
similarly hybridized (see Fig. 2). 

The spectrum in case (b) is determined by the follow- 
ing eigenvalue equation: 



lo{w 2 -4 cos 2 (fc/2) V 2 ) 



0. 



(12) 



so that there is a solution ui ~ 0. The appearance of such 
solution indicates that this configuration is energetically 
unstable. This conclusion is supported by the following 
calculations. 
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one gets for the other mean field solution:^ 
which has higher energy than the uniform mean field so- 
lution. 

Let us now consider fluctuations around the uniform 
mean-field solution. We write the g matrix in terms of 
two slow fields, as follows: 



g n = G{x) dl-P(x)+i(-l) n a-l(x) 



(17) 



where | l(x) \« 1 and G(x) belongs to SU(2). 

Neglecting fluctuations in V we get the following action 
on the lattice: 




(b) 

Fig. 2. Pictorial view of the hybridization of the chains. 
Cases (a) and (b) as described in text. 



The self-consistency equation determining V is: 

V = NJ < X i, m df m >, (13) 

which, upon introducing the one-dimensional tight- 
binding dispersion relation for the conduction bands, 
yields the following expressions for V, respectively, in 
the two cases: 



dk 



nj r /2 

2n J-*/2 Jsm 2 k + V 2 
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nj r /2 

~^ Jo 



dk 



(14) 
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One can see that V( a ) for the uniform mean field solution 
is larger than V(m ( by a factor of eih ), thus the ground 
state energy is likely to be lower in the uniform solution 
as we now proceed to show. 

The ground state energy is determined by the deriva- 
tive of the free energy with respect to the parameter V: 

OF V N 

gy = 7 + y (< Xi,m4 >m > + < d hmX t m >) , (15) 
which in the uniform case leads to: 



S = $+(d T + G+d T G)^ n + ^+e[id x + ad x G + G]$ r . 



+i\{x)(d 
+i(-l) n *+a-l(x)* n 



Xn,m 



(18) 



Let us now decompose the unhybridized fermion into 
slowly varying components, 



Vn,m=R(x)(i) n + L(x)(-i) n , 



(19) 



where x = na ( a is the lattice parameter), substitute 
this into the above effective action and integrate over 
fast fermionic fields x,<i with energies above the gap: ~ 
¥jj = Tk, and the field I. After the integration we get 
the following effective action for the flavour and the spin 
sectors: 



S 

S[PCF] 
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S[PCF] +S[GNl 



(20) 



^-Tr(d T G+d T G) + ^v 2 Tr(d x G+d x G), (21) 
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where v ~ vp and the spinor ip incorporates the R, L 
fields. Both effective actions have Tk as their ultraviolet 
cut-off. 

The action for G must vanish in the limit V — > 0. 
This will happen through the shrinking of the cutoff: 
Tk — * 0. 

The action in Eq.(^) corresponds to the Principal 
Chiral Field model on the SU(2) group. It is exactly 
solvableO; the excitations have a spectral gap 
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Fermionic model (P2|), describing the sain sector, is the 
0(2N)-symmetric Gross- Neveu modellla with the bare 
coupling constant g = n/16N. It-,is also exactly solv- 
able and all excitations have gapsli3. The mass scale is 
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~ T K e"H^2) _ T K e- 16 (24) 

This gap is numerically small despite the fact that it does 
not contain N. 

Thus, we see that the heavy fermion phase contains two 
modes, the spectra of which lie well inside of the single- 
particle gap. One mode is magnetic ( the one described 
by the Gross-Neveu model) and the other is not. 
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